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Abstract 

On the basis of pure gauge sector of Blaizot-Iancu equation, we derive kinetic equation 
of Boltzmann type, taking into account In + 2-colorless plasmon decay processes, n = 
1,2,... . Using so-called Tsytovich correspondence principle, a direct connection between 
matrix elements of the plasmon decay processes and certain effective current, generating 
these processes, is established. The procedure of calculation of matrix element for simplest 
four-plasmon decay is considered comprehensively. The limiting value of the plasmon 
occupation number (~ 1/g 2 , where g is a strong coupling) wherein all plasmon decays 
with n > 1, contribute to the right-hand side of the Boltzmann equation, is defined. The 
iterative method of calculation of matrix elements for higher decay processes (n > 1), is 
proposed, and a problem of their gauge-invariance is discussed. Proceed from the general 
reasons the problem of extension of suggested approach to the case of color plasmons, 
is considered. The explicit form of linearized Boltzmann equation for color plasmons is 
written out, and it is shown that this equation covariantly conserves the color current, 
resulting from color-plasmon number density. 
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1 Introduction 



This work consisting of three parts, is concerned with in-depth analysis of dynamics of 
boson excitations in hot QCD plasma at the soft momentum scale (~ gT, where T is a 
temperature of a system), within the hard thermal loop (HTL) effective theory. Although 
this range of momentum was studied in detail in the literature (see review of Blaizot 
and Iancu [|TJ), but nevertheless the dynamical processes of higher order in the coupling 
g connected with high- loop diagrams of an effective perturbation theory ||, remain to 
be practical analyzed. To such processes in particular we carry over the processes of 
the following type: the scattering of soft modes among themselves and the processes of 
scattering of two and more plasma waves by hard thermal particle. This and accompaying 
papers are concerned with research of the processes within a real-time formalism based on 
kinetic equation for soft excitations. For sufficiently high energy level of the soft plasma 
excitations (exact estimations will be given below) these higher processes of interaction 
play an important role identical to simplest scattering process including one soft collective 
quantum and hard particle, and is called nonlinear Landau damping. In the case of non- 
Abelian plasma, for boson excitations, the nonlinear Landau damping was studied in 
detail in Ref. ||. In it was shown that the nonlinear Landau damping rate is closely 
related to the damping rate obtained in Ref. Q in one-loop approximation of an effective 
perturbation theory. In the second part of our paper |J the extension of an approach 
developed in Ref. [[| to the case of scattering of arbitrary number of waves by hard thermal 
particle will be introduced, and exact estimations for oscillation amplitudes, wherein 
these higher processes are of the same order in g, will be given. The application of the 
theory developed there to the problem of energy loss of a color charged particle traversing 
the hot QCD plasma will be considered. In the third part the plasmon production by 
bremsstrahlung and processes of multiple scattering with emission of soft boson waves 
(longitudinal or transverse) will be studied, what can be useful for research of Landau- 
Pomeranchuk-Migdal effect in hot gauge plasma. 

Here, in Paper I we have restricted ourselves to study of the processes connected with 
self-interaction of the soft excitations. The study of these processes is to be preceded 
by consideration of more complicated scattering processes including soft modes and hard 
particles which contain self- interaction of the soft modes as a required constituent element. 

In our previous paper || a first attempt of research of dynamical process of nonlin- 
ear self-interaction of longitudinal colorless excitations (colorless plasmons) based on the 
Boltzmann equation for plasmon number density N l , was presented. The process of elas- 
tic scattering of two colorless plasmons for which the appropriate scattering probability 
was calculated, is considered. Unfortunately, method of derivation of collision term pro- 
posed there, is rather cumbersome and gives no way of its direct extension to calculation 
of collision term taking into account the scattering processes involved arbitrary (even) 
number of colorless plasmons. The present work is based on use of some different ap- 
proach allowing immediately to reduce the problem of calculation of the collision term 
to the problem of calculation of the scattering probabilities for appropriate scattering 
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processes. Above-mentioned assertion that for certain energy level of the soft excitations 
all processes of self-interaction of the plasma collective modes are of the same order in 
g, is stated more concrete in terms of the plasmon number density, and the equation 
of the Boltzmann type defines its space-time evolution. In present paper it shown that 
with the condition when the plasmon occupation number is of order 1/g 2 , the right-hand 
side of the Boltzmann equation represents a sum of an infinity number of terms having 
a sence of (2n + 2)-plasmon collision integrals, n = 1,2,.... The last circumstance is 
in reality a different restatement of well-known [|7|, |8| very nontrivial fact that for strong 
field A^{X) ~ T, all iV-gluon amplitudes with N > 2 contribute to the induced current 
j'^ nd [A] at the same order in the coupling. The HTL-amplitudes here represent coefficient 
functions in integrands for functional expansion of _7^ nd L4] in powers of gauge field A^. In 
this paper by introducing notions of interacting and free soft gauge fields, we somewhat 
deepen analysis of this circumstance. In particular, it shown that the induced current has 
the more physical content if it will be represented as an expansion not in interacting field 
Ap (as it actually takes place in Refs. [0, §), but in free field Affl. The coefficient functions 
of this expansion, called effective amplitudes, represent rather complicated combinations 
of HTL-amplitudes and an equilibrium propagator for a soft gluon. They preserve many 
features of usual HTL-amplitudes, but have more direct physical meaning - these func- 
tions set on mass-shell plasma excitations representing the scattering amplitudes of soft 
collective modes. 

At present there are a few methods of constraction of relativistic kinetic equations (see, 
e.g. review in Ref. 0). A more powerful and more convenient tool to derive relativistic 
Boltzmann equations from exact field Sc winger- Dyson equations, the so-called closed- 
time-path (CTP) formalism of nonequilibrium quantum-field theory, is considered. This 
general method is suitable to the same extent for deriving transport equations both hard 
modes and soft ones of plasma medium. From most papers close to the subject of our 
research, where for derivation of the Boltzmann equations for soft gluon longitudinal and 
transverse quasiparticles the CTP formalism is used, it should be mentioned the work by 
Niegawa ||. Our view on the problem of derivation of the Boltzmann equations for soft 
plasma modes is quite different from a view presented in Ref. ||. Since we are interested 
only with the processes of nonlinear interaction of the soft excitations among themselves 
and with hard termal particles, the use of quasiclassical approximation is justified. By 
virtue of this fact we consider that it has no need to use tools of nonequilibrium field 
theory which are general, but for our particular problems, are needlessly complicated. 

Our approach is based on the fundamental system of dynamical equations derived by 
Blaizot and Iancu ||, which is a local formulation of the HTL-equations of motion for 
soft fluctuating gluon and quark fields and their induced sources. These equations are 
obtained by means of a truncation of the Schwinger-Dyson hierarchy and already contain 
the maximum comprehensive information at leading order in g about dynamics of hot 
QCD plasma at soft momentum scaleQ Therefore this system is the most convenient and 

1 However it should be noted that the approximation used here may be unsuitable in attempt of 
calculation of effective amplitude of elastic plasmon-plasmon collision in the limit of vanishing momentum 
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optimal start point for derivation of transport equations for soft plasma excitations. In 
our previous paper || it was shown how starting from similar dynamic equations one can 
derive collision term for elastic scattering of two colorless plasmons. However, as was men- 
tioned above, the direct approach presented there, is very cumbersome and complicated 
for its extension to the processes of collisions with account of many colorless plasmons. 
For this reason in the present paper and accompanying ones, proceed from quasiclassical 
character of a problem, we initially assume that a structure of the collision terms are 
determined by Fermi's golden rule and thus we focus our efforts on calculation exclusively 
of the scattering probabilities, which in fact determine all dynamics of interaction of soft 
boson excitations. 

The following step is to extend analysis carried out in these papers on the fermion 
degree of freedom of plasma excitations in hot QCD plasma. Here, the series of new 
specific moments, which will be fully considered in the subsequent works, is arizen. 

The paper I is organized as follows. In section 2 preliminary comments, with regard 
to derivation of the Boltzmann equation, describing (2n + 2)-plasmon decay processes, 
are given. In section 3 all the conventions and the notations used in this paper, are sum- 
marized, and the nonlinear integral equation for gauge potential A^, playing key position 
in our subsequent research, is written out. Section 4 is devoted to explanation of the 
correspondence principal that enables one to establish a direct connection of the matrix 
elements of elastic collision of n + 1 colorless plasmons with certain effective current, 
generating these processes. Section 5 presents a detailed consideration of the Boltzmann 
equation for four-plasmon decay. In section 6 we estimate the typical value of plasmon 
occupation numbers, wherein one can rectricted the consideration to taking into account 
only contribution from four-plasmon decay or it should be considered all higher decay pro- 
cesses. In section 7 a complete algorithm of succesive calculation of the matrix elements 
defining (2n + 2)-plasmon decay processes in the temporal gauge, is presented. In section 

8 we discuss a problem of a gauge independence of these matrix elements. In section 

9 we give some speculation concerning extension of an approach proposed in this work 
to research of dynamics of colorless plasmons to the case of color plasmons. In closing 
section, when the main results were outlined, we briefly discuss the additional difficulties 
associated with appearance of the nonlinear frequency shift of longitudinal oscillations 
generated by plasmon collective interactions. 



2 Preliminares 

The main part of our work will be concerned with the constraction of an effecive kinetic 
theory for soft colorless longitudinal excitations, propagating in a purely gluonic plasma, 
with no quarks. In other words we assume that a localized number density of plasmons 

transfer, near to "forward scattering". 
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p,x) 



(K ab ) 



is diagonal in color space 



ab 



5 ab N l 



where a,b = 1,...,N^ — 1 for SU(N C ) gauge group, and we consider the change of 
the number density of the colorless plasmons N p as a result of their interactions among 
themselves. In section 9 only using some general reasonable assumptions, we made an 
attempt of an extension of the obtained results to the case of color excitations with non- 
diagonal plasmon number density. 



The dispersion relation for plasmons u l = u> l (p) = u l p is defined by 

Re*A _lz (t.;,p) = 0, 



(2.1) 



where 

*A- 1 W) 



r i + 



Hi 
p 2 



p 



X 



, F(x) = - 



In 



1 + x 



1 — x 



m9{\ - \x\) 



is an inverse resummed longitudinal propagator and = g 2 N c T 2 /9 is a plasma frequency 
square. From physical considerations it is clear that the notions of the plasmon and 
accordingly, plasmon number density have a meaning only with use of the condition 

oj 1 t 1 > 1, 

i.e. for the plasmon frequency that is much more than inverse lifetimes of plasmon 1/V. 
Generally speaking, this time is defined by just as the processes of plasmons scattering off 
each other, so by the scattering processes including soft modes and hard termal particles 
(hard thermal gluons in our case). For sufficiently large intensity of plasma excitations 
(exact estimations will be given in Sec. 6 of this paper and in 0), the processes of a first 
type can play the same important role in dynamics of a system as the second ones and 
even become dominant. Under this condition, plasma may be considered as involving two 
interacting subsystems: the subsystem of hard thermal gluons and the subsystem of soft 
plasmons, which exchange energy among themselves. As was mentioned in Introduction, 
here, we have restricted ourselves to detailed consideration of the processes of interactions 
in the plasmon subsystem. This fact is embodied in the structure of collision term in the 
plasmon kinetic equation proposed below. 

We expect the time-space evolution of scalar function N p to be described by 



dN l v 



V 



dNl 
<9x 



(2.2) 



where Vp = du l p /dp is a group velocity of the longitudinal oscillations. The generalized 
decay rate Td and inverse decay rate Tj in general case are (non-linear) functionals de- 
pendent on the plasmon number density. From here on such a functional dependence is 
denoted by argument of a function in square brackets. 
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We shall consider that the decay rate Td and regenerating rate T; can be formally 
represented in the form of functional expansion in powers of the plasmon number density 



oo oo 

UK\ = £ rf^livj], imjvJ] = £ if n+1) [iv p ], (2.3) 



n=l n=l 

where 



r (2n + l) [JV ,j = J dT (2n + l) ^ pn+i; ; ^ ^ (2 4) 



x(l + ^„ +1 )...(l + < n+1 ), 



r p"+i)[jVj] = / rfT( 2 ™ +1 ) «, 2n+2 (p, Pl , . . . , p„; p n+1 , . . . , p 2 „ +1 ) (1 + N l pi ) . . . (1 + JV*J 

(2.5) 



xJV ...M . 

Pn + 1 P2n + 1 



Here, w 2n+2 = w 2n+2 (p, Pi, • • • , Pn, Pn+i, • • • , P2n+i) is a scattering probability for process 
of elastic collision of n + 1 plasmons, and phase-space measure is 

/ dT^ = J (2vr)V 3 )( Pin - p out )S(E m - E Q Jf[^-, ( 2 -6) 



where 
and 



Pin = P + Pl + • • • + Pn, E in = J p + U) l pi + . . . + UJ p 

.1 , .1 



Pout — Pn+1 + • • ■ + P 2n+1, E ont — UJ pn+l + . . . + W p2n+1 - 

The Dirac ^-functions in Eq. ( |2.6D expresses the momentum and energy conservation in 
the collision process. In writing the expressions ( |2.4| ) and ( |2.5[ ), we have used the fact 
that the probabilities of direct and reverse processes are identical. 



We make some comments concerning Eqs. (|2.3|) - (|2.5|) . By virtue of the fact that pro- 
cesses of nonlinear interaction of odd number of the plasmons are kinematically forbidden 
by the conservation laws, in expansions of generalized rates ( [2.3|) we leave only odd terms 
in powers N p . Furthermore a structure of the terms of the expansions ( [2.4|) and (|2.5|) is 
chosen such that it taken into account the scatering processes only with equal number of 
the plasmons prior to interaction and upon it, i.e. the scattering processes of the following 

typefl 

g*+gi^g2 + S3> forn=l, 
g*+gi+g2~ ^Sl + Sl + St forn = 2, (2.7) 



The term the decay processes will be also used for the processes of ( |2.7| ) type 
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where g*, g*, . . . are the plasmon collective excitations. Such a choose of the structure of 
the decay and regenerating rates assumes that decay processes with unequal even number 
of incoming and outgoing soft external legs of the type 

g* ^ g*2 + gs + §4 + gs> forn = 2, 

etc, are kinematically suppressed as compared with "elastic" scattering processes fl2.7|) 
(speaking more exactly, kinematic regions of momentum variables in "elastic" and "in- 
elastic" scattering accessible by conservation laws, are not be covered each other, and 
a contribution of the last processes to the nonlinear plasmon dynamics to the order of 
interest is not important). 

The scattering probability u>2n+2 must satisfy the symmetry relations over permuta- 
tions of external momenta 

^2n+2(P, Pi, • • • , Pn] Pn+1, • • • , P2n+l) = ^2n+2(Pl; P, • • • , Pn', Pn+1, ■ ■ ■ , P2n+l) = • • • 

(2.8) 

= ^2n+2(Pn, Pi, • • • , P', Pn+1, • • • , P2n+l) = W 2n +2 (Pn+1 , • • • , P2n+l', P, Pi, • • • , Pn), 

which are a consequence of an indistinguishability of the colorless plasmons. The special 
consequences of symmetry properties ( [2.8| ) are the conservation laws of energy momentum 
and total number of colorless plasmons in the decay process (|2.7| ) 

£ = / (2^5 J ° N ° = mmt ' K = / (2^5 p N * = cmst ' 

3 Blaizot-Iancu equations. Correlation function of 
the boson excitations 

We adopt conventions of Blaizot and Iancu ||. We use the metric g^ u = diag(l, —1, —1, — 1), 
choose units such that c = ks = 1 and note X = (X ,X), p = (p = uj,p). The gauge 
field A u = A°^ a with — 1 Hermitian generators in the fundamental representation 
obeys the field equation 

[£>", F^{X)\ - ^ l n,n»A u {X) = Ja (X), (3.1) 

where D n = d a + igA n (X) is a covariant derivative, F av = F® u t a is a field strength tensor 
with F" = d^A" — d u A a ^ — gf abc A h ^A c v \ [, ] denotes a commutator and £ is a gauge 
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parameter fixing a temporal gauge. In what follows 4-vector n M will be identified with 
global 4-velocity of the plasma. 

The color current on the right-hand side of field equation ( |3.1| ) can be written as 

r dk 

j^X) = 2gN c j-— Vl ,5N(k,X), v = (1, v), v = k/|k|, (3.2) 

where 5N(k, X) = 5N a t a is a soft fluctuation in the gluon color density. We suppose that 
there is no external color current and/or mean color field in the system, and therefore 



an expectation value of the induced color current (|3.2j ) over the off-equilibrium ensemble 



equals zero, i.e. (j^{X)) = 0. The consequence of this fact is (8N(k,X)) = 0. 

On a space-time scale (gT)^ 1 the soft fluctuation 5N(k, X) satisfies the following 
equation: 

[v-D x , 5N(k, X)] = ig([v ■ A(X), 6N(k, X)))-g{v E(X) - (v • E(X))} (3.3) 

Here, E(X) = E a (X)t a is a chromoelectric field, E* = F i0 ; N(e k ) = l/(exp(e fc /T) - 1) is a 
boson occupation factor and e k = |k|. On the right-hand side of Eq. (|3.3|) we have entered 
the additional average terms, so that Eq. (|3~3| ) will be equal zero for the expectation value 
taken over the off-eqilibrium ensembles. The self-consistent system of Eqs. ( |3.1| ) - ( |3 .3| ) 
for a soft fluctuating field and their induced current presents a pure gauge sector of 
general system of self-consistent equations first obtained by Blaizot and Iancu in Ref. B. 



The equations ( |3.1|) - ( [3. 3D are solved by the approximation scheme method - the weak- 
field expansion. For this purpose first of all we expand the soft fluctuation of the gluon 
color density 

oo 

5N = J2 5Nis) > ( 3 - 4 ) 

s=l 

where index s shows that 5N^ is proportional to the sth power of A^. The expansion of 
a color current, corresponding to (|3.4j ) has the form 



00 r <7k 

U = E4 S) > ^ = 2 9N C J v,5N (s) . (3.5) 



Now we turn to the field equation (|3.1|), connecting the gluon soft field with induced 
color current j M . Let us rewrite this equation, explicitly separating free parts in (|3.1|) from 
interacting terms. We have 

d u {F„ v ) L - - J'i 1} = Jnl, - ig mA„, A v ] - (d v [Ap, A,])}- (3.6) 

-ig{[A v t (FM - ([A" 1 (F IV ) L ])}+ g 2 {[A\[A^A v ]] - ([A", [A,, A u ]])}. 

Here, indices L and NL denote linear and nonlinear parts of strength tensor and the color 
induced current with respect to A^. As in the case of dynamical equation (|3.3|) , we enter 
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the additional averaged terms on the right-hand side of Eq. ( |3.6| ) , which ensure vanishing 
left and right parts when an expectation value is taken. 

Substituting the expansion ( |3.4j ) into fl3.3| ) and collecting the terms of the same order 
in A^, we derive the system of equations for SN^ S \ s = 1, 2, ... , whose solution is easily 
obtained by Fourier transformation. Substituting such a derived expression we obtain the 
explicit form of the terms in a color current expansion ( |3.5|) 



jW(X) = t a jdpj^ a (p) e+ x , s = 1, 2, . . . , (3.7) 

where 

j«» = Tl^(p)A a »(p) } 

tf )a (p) = ^/^™ (p, - Pl , - P 2){A a ^( Pl )A a2fl2 (P2) - (A a ^( Pl )A a ^(p 2 ))} 

xSip-pi -p 2 )d Pl dp 2 , 

. . . , 

j ( ; )a (p) = j { g 3 - 1 J w™\::%(p, -pi, • • • , -Ps){A a ^( Pl )A a ^(p 2 ) . . . A a ^( Ps )- 

-A a ^( Pl )(A a ^(p 2 ) . . . A a ^( Ps )) - ... - (A a ^( Pl )A a ^(p 2 ) . . . A a ^(p s ))} 

s s 

x $(p-J2pi)Tl d P" 

i=l i=l 

Here, Ii^ v {p) is a soft-gluon self-energy and SF^" a * are usual HTL-functions |2|, |8[. 
Rewriting equation (|3.6| ) in momentum space and taking into account (|3.7|) we lead to the 
nonlinear integral equation for gauge potential A^, playing key position in our subsequent 
research 

oo 

•p-^CpKCp) = -r/M = - E J (s)a ^ A , -■■,*), (3-8) 

s=2 

where members of a series on the right-hand side of Eq. (|3.8| ) have a structure 

■•-,A) = j l 9 s - 1 ! t;;-/::;/.. (p, -pi, ... , - Ps ){A a ^ ( Pl )A a ^ ( P2 ) . . . A a ^ ( Ps ) 

(3.9) 

-A ai ^(p l ){A a2 ^(p 2 ) . ..A as ^(p s )) - ... - (A aiMl (pi)A a ' 2 ^ 2 (p 2 ) . ..A as ^(p s ))} 

s s 

xS(p-^2Pi)Hdpi. 

t=l i=l 

Here, the coefficient functions for s = 3 and s = 4 represent a sum of a bare 

vertices and corresponding HTL-corrections, and for s > 4 we rename *T^ S > = 
*V flu (p) is a medium modified (retarded) gluon propagator in a temporal gauge 

2 

*P^(p) = -P, u (p) *A t (p) - Q^ip) *A l (p) + ^-^—D^ip), (3.10) 
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where *A*''(p) = l/(p 2 - IP''(p)), U*(p) = \W u (p)P^ v (p), Ii\p) = IF* '{p)Q ^(p) . The 
Lorentz matrices in ( |3.10|) are defined by 



^V(p) = ^ ^ o—Qmp), QAp) = -21 \ — ' D MP) = 2~> 3 - n 

pZ U z {p) p z 



P 

u^ip) = j—ripv-u^ip-u)), u = p 2 u^-p ll {p-u). (3.12) 

Let us assume that we are in a rest frame of a heat bath, so that = (1, 0, 0, 0). 

At the end of this section we would like to introduce the correlation function of the 
soft-bosonic excitations, 

It W,P) = {Al\p')A b M). (3.13) 

The asterisk denotes the complex conjugate. The considered soft-gluon excitations are 
necessarily colorless by virtue of the fact that mean field (A°) or associated mean induced 
color current are assumed to be vanishing. Therefore, for the physical situation of interest, 
the off-equilibrium two-point function ( |3.13|) is diagonal in a color space. 

For the conditions of stationary and homogeneous of hot gluon plasma we have 

i;l(p , ,p) = 5 a % I/ (p')5(p'-p). (3.14) 

The off-equilibrium perturbations which are slowly varying in space and time lead to a 5- 
function broadering, and I^ v {p' ,p) depends on both arguments p and p'. Let us introduce 
Ifj,v{p'iP) — IfiuiPi Ap), Ap = p 1 — p with | Ap/p | <Cl, and insert the correlation function 
in the Wigner form 

I^(P, x) = J I^(p, Ap) e~ tAp - x dAp, (3.15) 

slowly depending on x = (£, x). In global equilibrium hot QCD plasma the oscillations of 
two types: the longitudinal and transverse ones can be extended. In this connection we 
define the Wigner function I^ u (p,x) in the form of an expansion 

I, u (p,x) = P^(p)P p + Q,ul l p , I™ = I^\p,x). (3.16) 



4 The correspondence principle 

The main purpose of this section is a derivation in an explicit form of the scattering 
probability W2 n +2 for process of elastic collision of n + 1 plasmons in hot gluon plasma. 
In our previous paper || we have derived the probability u>4 for a two-to-two scattering 
process of colorless plasmons. It was calculated by simple extracting all possible contri- 
butions to this process. However, use of such a direct approach in a general case of the 
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process of elastic scattering of n + 1 plasmons becomes ineffective because of awkwardness 



and complicated computations. In this case the method developed by Tsytovich |H], |1 1 
in the theory of nonlinear processes in electron-ion plasma, known as the correspondence 
principle, is more convenient for deriving explicit form of the scattering probability w 2n +2- 
We have already mentioned it in Ref. ||. Here, this approach will be developed as ap- 
plied to general case of (2n + 2)-plasmon decay. For non-Abelian plasma this approach 
is especially effective in the temporal gauge, when we have a closer correspondence with 
the electrodynamics of an ordinary plasma. The gist of this method is as follows. 

The change in the colorless plasmon numbers, caused by spontaneous processes of 
plasmon decays only, is 



/Fin 1 f)N l \ sp °° r 

+ af J / dTi2n+1) W2n+2{p > p " ; p " +1 ' • • • ' P2 " +l) < • • • ^ • 



This equation follows from (|2.2|) in the limit of a small intensity — > and using the 
fact that the occupation numbers iV* are larger than one, iV' + 1 ~ AT' In this case 
the change of energy of the longitudinal excitations is 



dS \ sp d / f dp 



dt I dt \J (2vr) 



\3 



E / ?^S / ^n +2 (p, Pi, ... , Pn, Pn + 1, P 2n + l)^ 1 . . . TV^ . 



71=1 ' 



(2 



On the other hand the value (d£/dt) sp represents the emitted radiant power of the longi- 
tudinal waves I 1 , which in turn is equal to the work done by the radiation field with the 
color current, creating it, per unit time 

r/2 

l l = lim — / /dxdt(E a (x,t)- J a (x,i)) = lim [dpduo (E a (p, u) ■ J a (p, uj)). 

t, v-^oo tV J J ' r, v^oo rV J 

-r/2 V 

(4.2) 

Here, V is a spatial volume of integration, E ai ('x,t) = — dA ai (x, t) /dt is a chromoelectric 
field in the temporal gauge, E ai (p, a>) and J ai (p, uj) are the Fourier components of field and 
current, correspondingly. Averaging over the time was made for elimination of oscillating 
terms in I 1 . For a transformation ( |4.2|) we use a relation 

r/2 



lim / e iujt dt = 2n5(uj). 

r^oo J 



-r/2 

The Fourier- component of a field E a (p, cj) = pE a (p, w)/|p| is associated with J a (p, uj) by 
the field equation 

^( P ,c) = ^*AW) (P- JQ (P^)) , (4.3) 
uj p 
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For weak-absorption medium, when Im*A 11 (p, u) — > 0, the propagator *A l (p,u) can be 
approximated in the following way 

Here, we consider that sign(Im*A _1 '(p)) = sign (a;). The symbol P denotes a principal 
value. Substituting Eq. ( |4.3|) into (|4.2j ) and taking into account reality of X 1 , we derive 

X 1 = -^lim^ ^ J dp W sign(w)<r'(p) (J£ a (p) j;(p))5(Re*A- 1 '(p)). (4.4) 

In derivation of the last expression, we present a longitudinal projector p ® p/p 2 in a 
Lorentz covariant form (Eqs. ( |3 . 1 1 ) and ( 3.12|) ) 



Since the nonlinear part of a current only is responsible for processes of plasmon decays, 
then it is necessary to set J a/t = J^ L in the expression ( |4.4|) , where J^ L is defined 
by equation (|3.8| ). We perceive the 5-function of the real part of inverse longitudinal 
propagator in the ordinary sense 

5(Re*A- 1 '(p, u)) = Z,(p) [6(w - «i) + % + wj,)], (4.6) 
where Zj(p) is the residue of the effective gluon propagator at the plasmon pole. 



Substituting the nonlinear current expansion Q3.8p into a correlation function in inte- 
grand of Eq. ( |4.4j ) we face the product of two series. However in this product it is necessary 
to leave only a sum of a product of terms having the same order in power of a potential 
A^. Then we have equal number of complex-conjugate potentials and the non-conjugate 
ones between the inside of the angular brackets of statistical averaging. The necessity of 
such a choose of terms was demonstrated by a direct calculation of the Boltzmann equa- 
tion for four-plasmon decay in Ref. 0. In this case only a desired 5-functions entering 
to integration measure expressing the energy and the momentum conservation of 

the four-plasmon decay, arise. We assume this for general rule, which is valid for decay 
process with arbitrary number of plasmons. Thus the emitted radiant power (|4.4f ), taking 
into account emission caused by plasmon decay processes only, can be presented in the 
form of expansion 

oo 

Z* = £z*«, (4.7) 

s=2 

where 

I Hs) = _ n ii m Mk W signHQ,>)(/| a >)^'( P ))J(Re*A- li (p)). (4.8) 

T, V — >oo TV J 
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In order to define the scattering probability W2 n +2, the correlation function on the 
right-hand side of Eq. (|4.8| ) has to contain terms of 2 (2n + 2)th order in a free field Aff' a . 
The required 2 (2n + 2)-order correlator yields the color current J^ a ( |3.9|) for s = 2n + 1, 
where the interacting field should be replaced by free field: A® — > However, here, it 

is also necessary to take into account the effects which arise from iteration of the currents 
j(s )a £ lower order, 2 < s' < 2n + 1. They give a contribution to the process of 
(2n + 2)-plasmon decay of the same order in coupling constant, as current J^ 2n+1 K The 
consideration of all similar contributions can be effectively presented as a replacement of 
a current J^ a [A] by certain effective current J^ a [A^]: 

J^ a (A,...,A)^J^ a (A^,...,A^) (4.9) 

= i g s - l J*tll\-;X^ -Pi' • • • ' -Ps){A (0)aiMl (Pi)A (0)a2M2 (P2) • • • A^ a ^( Ps ) 
-A^ a ^( Pl ){A^ a ^ 2 ( P2 ) . . . A (0)as ^(p s ))- . . . -{A {0)a ^ 1 (p 1 )A (0)a2IX2 (p 2 ) . ..A {0)asfls (p s ))} 

s s 
i=l i=l 

where the coefficient functions *r^ 1 i "^ 3 s (which in the subsequent discussion we shall call 
effective amplitudes as distinct from usual HTL-amplitudes *r^ 1 i "^ t 3 s ) represent highly 
nontrivial combinations of HTL-amplitudes and an equilibrium propagator for a soft 
gluon. Below an algorithm of their succesive definition will be presented. Here, the 
fact that a general structure of effective current Jjf' a is similar to a structure of initial 
expression (|3.9|) , is only important for us. 

Substituting ( |3.9| ) into Eq. ( |4.8| ) we obtain a sum of products of the correlators of 
different orders: 2sth-order correlator, product of pair correlators multiplied by (2s — 2)th- 
order correlators etc. The next step is the correlation decoupling of higher correlators 
in terms of pairs and next expressing (/l^ )^ )) in terms of plasmon number density iV'. 
However as was shown in Ref. the terms in the decomposition do not all associated with 
the processes of plasmon decays. The decomposition of averaging free field amplitudes 
into the correlators containing the pair of complex-conjugate potentials or one of the non- 
conjugate potentials between the inside of the angular brackets of statistical averaging, 
does not give a contribution to the processes of interest to us. The reason of this fact is 
similar to reason, by which we have dropped all crossed terms in a product of two series in 
initial expression ( f4.4|) . Here, it is not appeared the requred 5-functions in the phase-space 
measure, expressing the energy and the momentum conservation of the decay processes of 
plasmons. For this reason in substituting the expression (|4.9| ) into ( f4.8Q we can drop all 
the terms in braces in ( [4.9|) containing the averaging and lead to the following expression, 
instead of QE§) 

X 1 W = -vr^lim^ ^T^yr 2 9 2s - 2 J dpu S ign(u)8 aa ' Q^' (p) *f ]^(p, ~Px, . . . , -p.) 
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x *f^v; a ;,(p, -pi, . . . , -p;)(A*(°) ai ^(pi) • • • a*^^)^ ^^) . . . A^ a '^( P ' s )) 



x *W(p - $»^ 4 >(p ~ EK) *(He*A-"(p))n dpM. (4.10) 

In order not to overburden equations by symbol ' * ' we use sometimes a dagger ' f ' to 
denote complex conjugation. We write out decoupling of the 2sth-order correlator on the 
right-hand side of Eq. ( |4.10| ) in terms of pair correlators, which gives a contribution to 
the processes of interest to us. Suppressing color and Lorentz indices and employing a 
condensed notion, A\ = A^ ai (pi), we have 

(Af } . . . A*^A { ? ...A<S>) = { [(Af ] A { ?) . . . (Af^A^) + (perms, of 1', 2', . . . s') 

+ (perms. of 1, 2, ... s)} = s\{(Al i0) AP) . . . (A^Af) + (perms, of l', 2', . . . s')}. 

All terms on the right-hand side of the latter equality, that is obtaned from the first one by 
all possible permutations of arguments (p^, a[, % = 1, • • • , s, give the same contribution 
to the emitted radiant power I 1 ^ and therefore, here, we can write 

(A* (0) . ..A«°)AP . ..4?) = (s\) 2 (Af ] AP) . . . (A*® A®) 



i=l 

Here, in the last line we use a definition of the correlation function I^i(pi) in the equi- 
librium, Eq. (|3.14D , and leave only the longitudinal part in its expansion (|3.16|) . Now we 
substitute the obtained correlation decoupling into Eq. ( 4.10p and perform an integration 
over n!=i dp'^ 

2*W = -7r lim 47^"V d p/n^( 2 ^)1* (4) (P-Ep0l 2 wsign(a;)g^(p) (4.11) 



X 



n ^(ftrtefe -pi, • • • > -^)*rso>, -pi, . . . , -p.) n j*(ft)*(H**A-"(p)). 



i=l 



By a 5-functions squared here, we mean as usual \V1 

1 

Pi)J = 77 

i=l 



5 (4) (p-E^ 



(27T)' 



i=l 



To take into account weakly inhomogeneous and weakly non- stationary in the medium 
it is sufficient, within the accepted accuracy, to replace equilibrium spectral densities I l {pi) 
by off-equlibrium ones in the Wigner form ( |3.15| ): I l (pi) — > I l (pi,Xi), slowly depending 
on x. Futhermore, we take the functions I l (pi,Xi) in the form of the quasiparticle approx- 
imation 

I l (pi,Xi) = 45(^-0 + lL Pi 5(u i + u l Pi ), 4 = l\pi,Xi) (4.12) 
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and goes over from functions I p . to the plasmon number density 

iV^ = -(27r) 3 2^Zr 1 (p J )4. (4.13) 

Since external soft-gluon legs lie on the plasmon mass-shell, then account must be taken 
of the fact that conservation laws of energy and momentum (see, section 2) admit decay 
processes with even number of plasmons, i.e. it is necessary to set s — 2n+l, n — 1, 2, . . .. 

Taking into account above-mentioned, and Eqs. ( |3.12| ), ((OJ), ( |4.12| ) and ( |4.13| ) we 
obtain instead of Q4.11|) 

i „ r 2n+l „ J „2n+l j 2n+l 

x M"<"-'»"+'(p, -p„ . . . , -p 2 „+i)M«""<»»«(p, - Pl , . . . , -p 2n+1 ) (4.14) 

2n+l 

x {S^ - u l p ) + 6{u> + wj,)} II - O + ^ Pi %* + <)}■ 

1=1 

Here, we introduce 

M— 1( p,- Pl , ...,- P2n+1 ) =S 2n (^V /2 f-^U (4-15) 



it 2 (p) , 



x "Jf / ZKP!l ^ V2 f ^ *poai...03n+i {v _ pi _ p2 +1 ) | h „ 

representing the interaction matrix element for (2n + 2)-plasmon decay. 4-vectors 

on the right-hand side of (|4.15|) represent usual wave functions of longitudinal physical 

i ii 

gluon in temporal gauge, where factor Z z (p) provides renormalization of gluon wave 
function by thermal effects. 

Furthermore we multiply out terms in curly brackets in integrand of Eq. (|4.14 ). In 
terms containing the factor of N l _ p ,S(uJi + u l ) type we replace p« — ► — p« (in so doing we 
follow the rule u l _ p . — > —u l p ). By arguments presented in section 2 it should be left only 
the terms responsible for elastic plasmon scattering of (|2.7| ) type, and drop the remaining 
ones in obtained expression. In this case integrand in Eq. ( |4.14 ) has the form 

2(27r) 4 u;p{5 {4) (p+pi + . . .+p n -p n+ i- . . .-p 2 n+i)M* {a} (p, pi, . . . ,p„, -p n +i, • • • , -P2n+i) 

2n+l 

x M {a} (p, pi, . . . , p n , -p„+i, . . • , -p 2 n+i) + (perms, of signs)} ]J N l p .. (4.16) 

i=l 



15 



Here, for brevity we enter multi-index notation {a} = (a, a±, . . . , a2n+i) and fix a sign of a 
first argument of the matrix element . The factor 2 takes into account a contribution 
of the term S(u + u l p ). The permutation over sign means all possible placements of n + 1 
minus signs among external momenta pi, P2, . . . P2n+i- As the last step it should be 
performed a replacement of momenta such that arguments of 5-functions for all terms in 
braces of Eq. ( |4.16|) become equal to argument of 5-function entering into the first term. 
This enables us put outside the brackets of 5-function as a common factor. Confronting an 
obtained expression for x^ 2n+1 ) and corresponding term in expansion ( f4.1|) , one identifies 
the required probability W2 n +2- 

W 2 n+2(P, Pi, ... , Pn, Pn+l, • • • , P2n+l) (4.17) 
= M* {a} (p,pi, . . . ,p n , -p n+ l, . . . , -p 2 „+l)M {a} (p,pi, . . . ,p n , -p n+ l, . . . , -P2n+l) 
(Pi) (P»i)Pi 2 ) (Pl,-:P») 

+ M* {a} M {a} + Yl M* {a} M {a} + . . . + Yl M* {a} M {a} , 

l<il<n l<ii<i2<n 

where on the last line = M^(p, pi, . . . , p n , — p n +i, • • • , — P2n+i)- The summing 

(Pi) 

symbol denotes summing over all possible momentum interchange p^, 1 < %\ < n 

l<ii <n 

by momentum (— PjJ, n + 1 < ji < 2n + 1. The symbol Xa<\\<2<n analogously denotes 
a summing over all possible interchange of momenta pair (pi 15 Pi 2 ), 1 < %\ < i 2 < n by 
momenta pair (— pj i; — p j2 ), where n + 1 < j\ < j 2 < 2n + 1 etc. Let us consider as 
an example a case of four-plasmon decay, that corresponds to n — 1. By Eq. ( |4.17p the 
scattering probability has a form 

m(p, Pi; P2, Pa) = M*W(p, Pl , -p 2 , -p 3 )MW(p, pi, -p 2 , -p 3 ) (4.18) 
+ M*^(p, -p 2 , Pl , -p 3 )MW(p, -p 2 , pi, - P3 ) 
+ M*^(p, -p 3 , -p 2 , Pi)MW(p, -Pa, -P2, Pi)- 

Here, three terms imply the existence of three possible channels of this process, which 
change the plasmon number density iVp 

g* + gl # + gs, g* +g*2^g*l + gl, g*+g3-Sl+g2- 



The expression (|4.18 ) is suitable for checking on the symmetry conditions ( |2.8| ), which is 



imposed on the plasmon-plasmon scattering probability. 



5 The Boltzmann equation for four-colorless plasmon 

decays 

In this section, we review the main features of the Boltzmann equation for four-colorless 
plasmon decays which recently derived in Ref . . There are no new results to be reported 
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here, and they are given for coherent and complete subsequent discussion. By general 
formula ( 4.15|) the matrix element M aaiCt2a3 (p, pi, P2, P3) of four-plasmon decay is defined 
by the effective amplitude *^ a ^^ 3 (p, Pi, P2, P3) that in turn is a coefficient function in 
integrand of an effective current J^^fA^ ^] (Eq. Q4.9| )). Such a problem is reduced to a 
definition of explicit form of Jj^'' 1 ^ ']. For this purpose we use a basic field equation 
( p.8| ), where on the right-hand side we leave the terms up to the third order in powers of 

A. 

*V- lfiV (p)A a v (p) = -J (2)afl (A, A) - J® a **(A, A, A), (5.1) 

where, 

J^ a {A,A) = ^/T™(p, - Pli - P2 ){A a ^( Pl )A a ^(p 2 ) - (A a ^( Pl )A a ^(p 2 )))} 

x 5(p -pi - p 2 )dpidp 2 , (5.2) 

and 

J^ a (A, A,A) = ± g 2 J *T a ;^ 3 (p, - Pl , -p 2 , -p 3 ){A a ^{ Pl )A a ^{p 2 )A a ^(p 3 ) 

-A a ^{p l ){A a ^{p 2 )A a ^{p 3 ))~{A a ^{p l )A a ^{p 2 )A^ 
The function 

*r™ (p, -pi, -P2) = -*f aaia2 T W1W ( P , - Pl , - P2 ), 
*CX 3 afe -pi> -p*> -»») = -r ib f ba2aA *r» m ( P , - Pl , - P2 , -ps) 

_ f aa 2bfbai a s T^^fo -ft, ~ Pl , ~Ps), 

are three and four-gluon HTL-amplitudes accordingly. 

The nonlinear integral equation ( |5.1| ) is solved by the approximation scheme method. 
Discarding the nonlinear terms in A on the right-hand side of Eq. ( |5.1|) , we obtain in the 
first approximation 

*v- 1 ^(p)a:( p ) = o. 

The solution of this equation, which we denote by A^ a (p), is the solution for a free field. 

Furthermore keeping the term, quadratic in field on the right-hand side of Eq. Q5.1| ), 
we derive 

*f>- ltiV {p)A a v {p) = -J (2) ^(A (0) ,A (0) ), 

where on the right-hand side we substitute free fields instead of interacting ones. The 
general solution of the last equation is given in the form 

A;(p) = Af\p) - *V, u {p)j^{A^\A^), 
where J^{A^\A^) = J^(A^,A^). 
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The following term in the expansion of the interacting field is defined from equation 
*f)- 1 ^(p)A a u (p) = -J^{-*VJ^ 2 \A^\A^),A^) (5.3) 

-jWw(A®\-*VJV\A®\A<®)) - J^(A(°\A(°\A^). 

The first two terms on the right-hand side of this equation represent iteration of the 
currents Jjf ^ a of lower order (in this case, s 1 = 2) mentioned in previous section, whose 
contribution is of the same order in the coupling g as J^ a ^ (A^°\ A^°\ A^). Using explicit 
expressions for currents (|5.2|) , after cumbersome algebraic transformations, we obtain 
the form of interacting field from the Eq. ( |5.3| ) with accuracy required for our further 
calculations 

A*(p) = Af a (p) - *V lliy (p)J^ au (A^°\A^) - *V flu (p)J^ au (A^,A^°\A^). 
Here, third-order color current on the right-hand side is defined by expression 

^3)^(0)^(0)^(0)) 

= -P3){^ ( ° )aim (Pl)^ (0)a2/12 (P2)A(°) a3 ^(p3) (5.4) 

_^(0)a lW ^ ^(0)a 2M2 ( p2 )^(0)a 3 M 3 _ ^(0)a lMl ( pi )^(0)a 2 ^ 2 (p 2 ) A^ a ^ (p 3 ))} 

x5(p-p 1 -p 2 - p 3 )dp 1 dp 2 dp 3 , 
where the effective four-gluon amplitude *r^^ 3 (= *f( 4 )) is defined by expression 

*f™%(P> - p i> ^ ^ = -f aaib f a2a3 ^^(p, - Pl , -p 2 , -p 3 ) (5.5) 

- t a2b f baias *f ^(p, -pa, -p 1; -p 3 ). 

The color factors in Eq. ( |5.5|) are multiplied by purely kinematical coefficients, which we 
shall call partial effective amplitudes or effective subamplitudes^, and defined as follows 

^wiwsiP' ~Pi> -Pii -Ps) = * r /w2/«(P; ~Pu ~P2, -Pa) 

- *r MMll ,(p, -pi, -p + Pi) *V uu \p2 + Pz) *r^ 2A13 (p 2 + p 3 , -P2, -Pa) (5-6) 

- *r w3 „(p, -pa, -p + p 3 ) *D w '(pi + p 2 ) T^ M2/11 (pi + p 2 , -p 2 , -pi). 

For complete picture we leave the term (A^A^A^) in braces under the integral sign on 
the right-hand side of Eq. (|5.4j ). It is clear that this term is equal to zero, because A^ 
represents the amplitudes of entirely uncorrelated waves. The expression ( |5.4j ) represents 

3 Here, we use the terminology accepted in theory of multi-parton hard processes (see, e.g. Q), where 
the expansion of (5.5) type in calculation of the tree QCD amplitudes in the high-energy limit, are also 
arizen. 
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desired effective current J^ a [A^], whose coefficient function by a formula ( 4.15 ) defines 
the matrix element for four-plasmon decay 




M a ™(p,- Pl) -p 2 ,-p3) =g 2 [ ) ( j (5.7) 




v?{pi 



on— shell 



The obtained matrix element has a simple diagrammatic representation drawn in Fig. 1. 
The black square here, denotes an effective amplitude *r( 4 ). The first term on the 
right-hand side Fig. 1 represents a direct interaction of four plasmons, that induced by 
usual four-gluon HTL-amplitude. The remaining terms are connected with plasmons 
interaction, induced by three-gluon HTL-amplitudes with intermediate virtual oscilla- 
tions representing s- and t-channel contributions respectively (effective subamplitude 
*r m2AtlM3 Go, — P2, — Pi, — P3) m Eq. (|5.5|) contains also ^-channel contribution, which we 
have not drawn in Fig. 1). 



With the matrix element ( |5.7| ) in hand it is not difficult to write the Boltzmann 
equation describing the plasmon-plasmon scattering. Here, we have restricted our con- 
sideration to linearized version of this equation. For this purpose we assume that the 
off-equilibrium fluctuation is perturbative small and write the number density of colorless 
plasmons as 

N l p = N l eq (p) + 5N l p , 

where N l eq (p) = (e^p/ T * — is the Planck distribution function and T* is a certain 
constant, which can be interpreted as a plasmon gas temperature in the statistical equi- 
librium state. Parametrizating off-equilibrium fluctuation of the occupation number 5N p 
as follows 

dN l (v>) 1 

*K = ^Wj = ^<(p)«(p) + l)W p , (5.8) 



p 



we derive from Eqs. ( |2 . 2|) - ( fE5|) for n — 1, and Eqs. (|4.18| ) and ( |5.7|) , after simple color 



algebra, a linearized Boltzmann equation for function W p 



dt p <9x J (2tt) 3 (2tt) 3 (2tt) 3 



P2 P3' 



x , ^ N <( Pl )«(p 2 ) + l)«(p 3 ) + l) ,, Q , 

x 5(p + Pl - P2 - P3 ) mp) + 1) ^ 

xw(p, Pl ; p 2 , p 3 ){W l p - W p2 + W pi - W l p J, 

where the function 

w(p,Pi;P2,P3) = 3# 4 iV c 2 ( - \ ) II ( o 7 K -2f \ ) {l*r(p ; -P2,Pi,-P3^ 12 
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+ I *f (P,P1, ~P3, -P2)\ 2 + Re(*f(p, -p 2 ,Pl, -P3) *f t (P,Pl, ~P3, -P2))} Ion-shell (5.10) 

is the probability of four-colorless plasmon decay. Here, we denote 

*f{p,p 1 ,p 2 ,p 3 ) = *f w ^^ 3 (p,pi,p 2 ,p 3 )^(p)u w (pi)^ 2 (p 2 )£t M3 (p3). (5.11) 



The structure of the expression on the right-hand side of Eq. ( |5.9|) is just as the structure 
of the expression for the case of the Boltzmann equation for hard gluons [Q. However, 
here, we have more complicated and cumbersome expression for scattering probability of 
two soft quasiparticles involving the resummed vertices. Besides, the momentum transfer 
q = p — P2 = P3 — pi is of the same order as the momentum of the soft quasiparticles. In 
the last case it is impossible to expand functions W p2=p _ q and W p3=pi+q as power series 
in momentum transfer q, as it takes place for hard gluons [13]. 



Using the linearized Boltzmann equation (5^) it is not difficult to estimate the order of 



the lifetimes of plasmon t 1 caused by the process of four-plasmon decay. Considering the 
plasmon gas in thermal equilibrium with hard particles from the heat bath, i.e. T* ~ T, 
and using Eqs. ( |5.9|) and ( pMOj ), we obtain 

i ~ g*N 2 c T. 

Thus necessary condition uj 1 t 1 ^> 1 (section 2), such that excited state of plasma can be 
described in terms of the plasmon occupation number iV p for g -C 1, is high accurate, 
since here, we have uj r ~ l/g 2 - 



At the end of this section we also represent some properties of the effective subam- 
plitude *r wn/J2At3 (p, — pi, — p 2 , — P3) defined by Eq. (|5T6|). In particular, deriving Eqs. ( p.9| ) 
and ( |5.10| ) we have used two relations which satisfy an effective subamplitude 

*f WtlM2At3 (p, "Pi) -P2, -P3)+*rV 2mM3 (p, -P2, "Pi, -P3)+*f WlAt3 M2(P ; -Pit -P3, -P2) = 0, 

(5.12) 

*f MWM2M3 (p, -Pl, -p 2 , -P 3 ) = *f W3M2Ml (p, -p 3 , -P2, -Pi). (5.13) 

Their correctness may be verified by direct calculation using known properties of HTL- 
amplitudes f2j, entering into the definition of *T ( - 4: \ These relations show that the effec- 
tive amplitude *n 4 ) possesses some properties of usual four-gluon HTL-amplitude. As 
we shall show below this statement takes place for the effective amplitudes *T ( - S ^ for any 
values s > 4 (the case when s = 3 is trivial by virtue of *r( 3 ) = T^). Here, the exception, 
for example, is the Ward identity, which for *r( 4 ) has not a standard form. 



6 Characteristic amplitudes of the soft gluon field 

In this section we shall estimate for which typical amplitude of the soft gluon field, the 
contribution of four-plasmon decay process to a generalized decay and regeneration rates, 
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TdfiVp] and I\[iV p ], will be leading and for which value all terms in the expansions ( |2.3| ) 
will be of the same order in g. At the last case two problem arise: the computation of 
explicit form of matrix elements of decay processes of high orders, for n > 1, and the 
problem of summing ( |2.3|) . Recall that in this work we consider the collective processes 
only at the soft momentum scale for u, |p| ~ gT. 

First of all we estimate an order of matrix element M aai '" a2n+1 . We use general ex- 
pression fl4.15|) connecting matrix element for (2n + 2)-plasmon decay with the effective 



amplitude *f ( 2n + 2 ) i n which the kinematical factors corresponding to external soft-gluon 
legs are single out. In the soft region of the momentum scale the following estimation 
results from this expression 



I2n + 1 
MMl ■■■A t 2n+1 ' 



Here, the symbol ~ denotes that the order of value from the left in coupling constant 
equals to the order of value from the right. The order of effective amplitude *r^ 1 1 '''" 2 2 ™ +1 1 
can be estimated from arbitrary tree diagram with amputate 2n + 2 soft external legs. Let 
us consider, for example, the diagram drawn in Fig. 2. Four-gluon vertices are of order 
g 2 , and HTL-resummed propagators are of order l/(gT) 2 . From simple power counting 
of the diagram it follows an estimation 

*-paai...a2n+i ^ 2nj\rn ^ 

and thus 

M oai...a 2 „+i ^ a 2n N n 1 (Q I) 

Furthermore integration measure has an estimation 

dT (2n+l) _ (gTf n -\ (6.2) 

Let us define an order in g of the plasmon density. For this purpose we use initial relations 
connecting the plasmon number density iV p with two-point correlation function ( ^.13[ ) of 
soft-gluon field. Here we have 



iV p = -(27r) 3 2 Wp Zr 1 (p)/ i 



i 

p- 



where in turns the spectral density I l p is defined from relation I l p = I l p 5{uj — u l p ) + I 1 _ v 5{oj + 
LJ l p ), and finally function I l p is defined by two-point correlation function 

(A^p')A u (p)) ~ Q^l^p' -p). 

Taking into account above-mentioned, we derive the following expression for estimation 
of the plasmon number density 

N l p ~ u l p Zi\p) [S( U - u l p )} ^ [S(jpf - p)} ~\a»{p>)A u (j>))- (6-3) 

Let us consider more typical two values of amplitude of a gauge field A^ at the soft 
momentum scale . Using Eq. ( |6.3D we have: 
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1. let \A,{X)\ ~ T (\A,(p)\ ~ 1M<?T) 3 ), then iV p ~ 1, (6.4) 

2. let ~ ^ T (I4*(p)I ~ UV9(gT) 3 ), then iV p ~ i. (6.5) 

For latter value of amplitude of the gauge field, when j<1, the occupation number 
Np is large, and the use of classical description is justified. The value of the amplitude 
of a gauge field (~ \f§T) is more probably at the momentum and energy scale of the 
plasmon modeQ. Power counting of the decay and regenerating rates, ( |2.4|) and fl2.5|) with 
regard to ( |6.1| ) and ( |6.2|) gives the following estimation 

r (2n+l) aa> _ g^ gT )N^(N l p ) 2n+1 . 

If we now set iV p ~ — , a>0, then from the last expression it follows 

r (2n+l) aa' _ g 2n{2-a) g X-a N 2nj, _ ( g g) 

For small value of oscillation amplitude (6.5) we have 

/ r (2n+l)aa'\ ^ g^^nj, 

From this estimation it can be seen that each subsequent term in the functional expan- 
sions (|2.3| ) is suppressed by more power of g 2 and here, we can only restrict ourselves to 
first leading term, describing four-plasmon decay process. Besides, here, the use of the 
linearized Boltzmann equation for colorless plasmons, i.e. Eq. (|5.9|) is justified , such the 
Planck distribution, relative of which the departure of the plasmon number density is 
defined, is of order 

<(P) ~ — ~ -• 
9 

In this case it is said that a theory of a plasmon-plasmon interaction for small amplitude 
of soft excitations is linear (linear amplitude regime), and nonlinear effects connected with 
nonlinear terms over off-equilibrium fluctuations of plasmon occupation number <5iV p can 
be treated as perturbations. 

A situation is qualitatively changed, when a system is highly excited. In a limiting 
case of a strong field, ~ T, for a = 2 from the estimation (6.6) it follows that 
functions p^ +1 ) aai j s no ^ ^ depend on n. All terms in the expansions ( |2.3| ) become of 
the same order in magnitude, and the problem of resummation of all relevant contributions 

4 Blaizot and Iancu [yj showed in the special case when the soft fields were thermal fluctuations at soft 
scale gT, that their typical amplitudes would be of the order |A^(X)| ~ \f§T. A similar estimation of a 
value of the classical fields at the soft scale was shown by Nauta and van Weert in Ref. [^5j based upon 
a somewhat different context. 
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arises. It is evident that a procedure of a linearization of a kinetic equation for plasmon 
occupation number in this case becomes unacceptable and here, we lead to truly 
nonlinear interaction theory of the soft excitations in hot QCD plasma. 

It is noted an interesting analogy with phenomenon of so-called parton saturation, 
which is of great interest in the physics of nuclear and hadronic processes in the regime 
where Bjorken's x becomes very small lfL6| , |17|. Saturation is expected at gluon phase-space 
density of order 1/g 2 (this is just maximum value of the plasmon occupation number that 
in our case corresponds to maximum value of an oscillation amplitude of soft gluon fields: 

~ T). Using this fact McLerran and Venugopalan |18| suggested a classical effective 



theory to describe the gluon distribution in large nuclei and valid for some range of x. The 
consequence of a development of this theory is a construction by Jalilian-Marian, Kovner, 
Leonidov and Weigert of renormalization group equation JOJ (the JKLW equation) that 
predicts the evolution of the gluon distribution function. In the saturation regime, where 
A^ ~ 1/g, a contributions of all orders in the strong background fields which represent 
the condensate, should be taken into account. In the low density, or weak-field limit, this 
equation is linearized and reduced to known equations. In our case the role of the JKLW 
equation fulfils a kinetic equation (|2.2|) that is complete nonlinear in a limiting case of a 
highly excitations of a soft-gluon field A^ ~ T and is reduced to a linearized Boltzmann 
equation (|5.9|) in the weak- field limit. 



7 The matrix elements for (2n + 2)-plasmon decays 

As shown in section 4, a calculation of matrix elements M aai "' a2n+1 which are responsible 
for processes of (2n + 2)-plasmon decays reduces to calculation of certain effective currents 
J^ 2n+1 ' ) being functionals of 2n + 2 free soft-gluon fields A^\ n = 1, 2 . . .. The coefficient 
functions in integrand of effective currents ( |4.9|) by equation ( |4.15 ) define matrix elements 



of decay processes in tree approximation, which allowed by conservation laws of energy 
and momentum. In section 5 an example of construction of such an effective current 
in simplest case of four-plasmon decay was presented. However as we have shown in 
previous section, it becomes necessary to take into account the contributions of higher 
decay processes, i.e. subsequent terms in expansions ( |2.3| ) when the energy level of the soft 
plasma excitations increases. Here, we consider the problem of determination of explicit 
form of matrix elements for higher decay processes leaving aside more subtle questions 
connected with summing and convergence of a series (|2.3|). 



As we have seen in section 5 the effective currents appear in the solution of the non- 
linear field equation (|3.8| ), that defines interacting soft-gluon field A^ in the form of a 
functional expansion in free field Ajp . Here, we shall extend this approach for determina- 
tion of the effective amplitudes for s arbitrary values. Let us rewrite here, the field 
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equation ( |3.8| ) for convenience of further references 



•^(pK(p) = = - E j(s)aM (A • • • (7.i) 



s=2 

where 



jW fl (A . . . , A) = ^g- 1 f m T%;;£(p, -pi,..., - Ps )A a ^( Pl )A a ^(p 2 ) . . . A-"- (p.) 

(7.2) 



x<*(p - n ^ 

i=l i=l 



On the right-hand side of Eq. ( [7.2|) we discard all terms, containing the thermal av- 
eraging (see discussion following by Eq. ( J4.9| ) in section 4). The nonlinear field equation 
( [7.1D can be formally solved perturbatively (at least in the weak-field limit) order by order 
in the coupling constant. We expand 

oo 

^(p) = E4 s_1)fl (p). (73) 

s=l 

where A^ 1 ' 01 ^) is a contribution of order g s ~ l to the soft-gluon field. For s = 1 we 
have a solution for a free field. Substituting an expansion ( |7.3| ) into Eq. ( |7. 1| ) we obtain 
iterative solutions of higher order in s: 

4 1)a W = -V,u( P )J {2)au (A(°\A^), 

Af a { P ) = -*V^{p)J® w (A®\A®,A(% 
Af a (p) = -^ F (p)J (4)a "(A (0 U (0 U (0) , ^ (0) ), (7.4) 

4 s " 1)a (p) = -*^(p)J ( ^(A(°),...,A(°)), 



Here, the effective current J^ av is defined by Eqs. ( |5.4|) - ( |5.6|) , the current J^ au is defined 
by next iteration 

j(A)au (4(0) ^(0) ^W^W) = J^Wy^, A (0) , A (0) , A (0) ) 
_ I j(3)a, ( _*pj(2) (A (0) j A (0) ); A (0) ; ^(0)) + j(3)a, (A (0) ) _*p J(2) (A (0) ^ ^(0)) ^ ^(0)) 

+ J( 3 ^(AW,AW,-*PJ( 2 )(A( ),^))} (7.5) 

_ | j(2V(_*pj(3)( A (0) j A (0 }) A (0) ); ^(0)) + j(2)a, (A (0) ; _*p J(3) (^(0) ^ ^(0) ^ ^(0))) j 
_j(2)a, ( _*pj(2) (A (0) )A (0) ); _*pj(2) (A (0) )A (0) )) 
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etc. A similar iterative procedure of computing perturbative solutions of the classical 
Yang-Mills equation in various modifications is often used in solving specific problems. 
From more closely works to subject of our research mention may be made of the works of 
Kovchegov and Rischke [20] , Matinyan, Miiller and Rischke concerned with the prob- 



lem of the classical gluon radiation in ultrarelativistic nucleus-nucleus collisions. Further- 
more, in Ref. p2] Bodeker used iterative solutions of the classical Yang-Mills equations 



for soft modes of the gluon field in deriving the Boltzmann equation for hard trans- 
verse gluons. Finally in the construction of the correlation functions for ultrasoft fields 

fa;, 



|p| ~ g 2 T), a similar iterative procedure was used by Guerin in Ref. |23 



However such a direct approach for determination of an explicit form of the higher 
effective amplitudes J T (s ), s>3, becomes very complicated and as consequence ineffective. 
In particular this associated with a necessity of additional symmetrization of integrands, 
which becomes extremely intricate with increasing order s. It is easily seen, for example, 
in attempt to determine the effective amplitude *r^ 5 ^ immediately from the right-hand 
side of Eq. (|7.5|). Here we proposed a somewhat different approach of calculation of *T (s ). 



It is not transparent as the approach based on iterations (|7.4|), but it allows us to avoid 



many intermediate operations and somewhat automate calculation procedure. 

The calculating algorithm is based on a simple idea. The nonlinear current J^[v4] 
has two representations: by means of free and interacting fields, which must be equal each 
other 

oo oo 

JmM = E J WaM (A (0) , . . . , A (0) ) = £ J {s)a>M (A, . . . , A). (7.6) 

s=2 s=2 

Here, the interacting fields on the right-hand side of the last equation are defined by 
expansion 

oo 

A a »(p) =A^ afl (p)-*V^'(p)J2J { J )a (^ \...,A^), (7.7) 



s=2 



where 



j^ a (A^, . . . , = - <r 1 J t;;-:::;,; (p, -pu ■ ■ ■ , - Ps )a^ ( Pl ) . . . a®-- ( Ps ) 

(7.8) 

s s 
i=l i=l 

Substitution of Eqs. ( |7.7| ) and (|7.8j ) into (|?T6|) turns this equation into identity. Now we 



functionally differentiate left- and right-hand sides of equality ( f7.6[) with respect to free 
field A^ afl considering Eq. ( |7.8j ) for differentiation on the left-hand side and Eqs. ( [7.2| ), 
( [7.71) and (|7.8| ) for differentiation on the right-hand side, and set A^ afl = after all 
calculations. The required effective amplitudes will appear both on the left-hand 
side and on the right-hand side. However main effect is that the effective amplitudes on 
the right-hand side will have at least one less external soft leg than on the left-hand side. 
This enables us to calculate them in a recurrent way. Below we shall give a few examples. 
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The second differentiation of a nonlinear current yields 



= *rz%(p> "Pi. -P2)^\p -Pi- Pa) 

A(°)=0 



5A(°)°i'*i(pi)<JA( ) a2 ^ 2 (p : 

= *KZ7M-Ph-P2)s (A) (p-Pi-P2). 



This example shows that the effective amplitude of three-plasmon decay coincide with 
usual three-gluon HTL-amplitude. This decay process is kinematically forbidden and 
hence the (^-function no support on the plasmon mass-shell. 

The first interesting case arises in calculation of the next derivative. It defines the 
effective amplitude for four-plasmon decay 



A(°)=0 



SAW ai/n (j,j)M(0) (p 2 )SAM a ^ {p. 

= *r7JXZM-Pi>-Pz>-Ps) 5i4) (p-Pi-Pz-Pz) = [*r^X 3 3 fe-ft>-P2,-P3) 

-{r?; 2 (P, -P+P2, -P2) *^>i+P3)*f^ 3 (pi+p 3 , -p 1; -p 3 ) + (circular perms, of 1, 2, 3)} 

X5 (4) (p-p x -p 2 -p 3 ). 

By using Jacobi identity for anti-symmetric structure constants, the right-hand side of the 
last expression can be led to the expression ( |5.5| ) obtained before. It can be viewed as an 
expansion of an effective amplitude *r""*" 2 2 "f 3 in terms of two independent combinations 
of a product of two structure constants 

jaa^b jbazaz jaaib jba\as (7 9) 

as it occurs for usual four-gluon HTL-amplitude || . A similar expansions take place also 
for higher functions. Let us consider the fourth derivative of a nonlinear current 

with respect to A^ a ^, defining a five-point effective amplitude. After rather cumbersome 
but direct calculations we obtain 



s 4 JnlM(p) 



A(°>=0 



<L4(°) (pi) . . . <L4(°) a ^ (p 4 ) 

= *^:7M -pi. • • • . -p 4 )5 (4) (p - Eft) = [*C££(p« -ft. • • • > -ft) 

i=i 

- { r ™(ft -ft - ft. -ft. -ft) *^ w '(ft + ft)* f I'ZM + ft. -ft. -ft) 

+ r^ 3 (P. "Pi, -P2 - P4, -P 3 ) *^'(p 2 + P4)*f5g^(P2 + P4, -ft, -P4) 

+ (circular perms, of 2, 3, 4) | 
-{ r ?A (P. -P + ft- -Pi) (pa + p 3 + p 4 )T S^Su(p - Pi, -p 2 , -p 3 , -P4) 



(7.10) 
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1 



+ (circular perms, of 1, 2, 3, 4) j 

+{r f; x ( P , - Pl - p 3 , -Pa - pa) *v vv \pi + P3)*f b ^M + p». -pi> -ft) 

x *P AA '(p 2 +p 4 )*f™ 4 (p 2 +p 4 , -p 2 , -p 4 ) + (circular perms, of 2,3,4) }]<5 (4) (p - 

i=i 

The obtained expression as it stands is very complicated and has involved color structure. 
By direct algebraic transformations, using the Jacobi identity for structure constants and 
relations (gig) and ( gig) it can be represented in the form of an expansion in terms 
of six independent combinations of a product of three structure constants, which we are 
choosen as follows 

j?aaibi jb\a,2b2 y&2i3«4 jaa\b\ jb\azb2 ^20,20,4, jaazb\ ybia^fe yfeaiiM 
j?aa,2bi j?bia,3b2 yfeaiey jaasbi jb\a\b2 ^620204 jaa2b\ jb\a\b2 jb2a-$a.A, 

Here, as in basis ( [7.9| ), we fixed the first and the last color indices. The expansion *r( 5 ) 
in terms of a basis ( [7.11]) reads 

*r^v: a ;> -pi, ■ -pa) = r^f^f^ *fv 1 ... w (p, - Pl , . . . , -pa) (7.12) 

+ (perms, of 1, 2, 3), 

where a sum is over all 3! permutations of the external gluons 1, 2 and 3, and effective 
subamplitude is defined as 

*f Mm ...^(p, -p u . . . , - P a) = *r Wil ... M4 (p, -pi, . . . , -p 4 ) 

- { r wl /^(P, -Pi, ~P2, ~Pa ~ Pa) *V u "'(p 3 + p 4 )*f „ W4 (p 3 + Pa, ~Pa, ~Pa) 

+ (circular perms, of 2, 3, 4) | 
- {r wi ^(p, -pi, -p + pi) *P w '(p - pi)*f v > milzi xSj> - Pi, -p 2 , ~Pa, ~Pa) (7.13) 
+ r M ^ 4 (p, -p + p 4 , -p 4 ) *P^'(p - p 4 )*f u'^^ip - Pa, -Pi, ~P2, -pa)} 
+ r Ml/A (p, -Pi - P2, -Ps - Pa) *V uu '(p l + p 2 )*t y ' M1M2 (pi + p 2 , -pi, -p 2 ) 
x *V xx '(p 3 +p 4 )*f A' M3M4 (P3 +P4, -Pa, -p 4 )- 

The examples considered above suggest that a color structure of the effective am- 
plitudes entirely coincides with color structure of usual s-gluon HTL-amplitudesf], 
derived by Braaten and Pisarski [23] for an arbitrary number of external soft-gluon legs. 

5 It is clear that if this statement is really true (for the present we cannot rigorously prove this fact) 
then the algorithm of calculation of is not optimal. It takes its additional modification, which 

enables one to derive at once the effective subamplitudes instead of the total effective amplitudes. 



27 



Besides, by direct calculation it can be shown that an effective subamplitude ( |7.13| ) sat- 
isfies two relations which hold for usual five-gluon HTL-amplitude, namely 



* T ^W2^(P, -Pi, -P2, ~P3, -Pi) + * r ww^2M4(P> ~P3, -Pi, ~Pi, ~Pa) (7-14) 

+ 'fwiiw/ulP) ~Pi, ~P3, -Pi, ~Pi) + *f WlM2Al4M3 (p, -Pi, -ft, ~P4, -Pa) = 0, 

*f' w * 1AtaWA44 (p, "Pi, -P2, -P3, -P4) = - *f M(tMM3M2Atl (p, -p 4 , -P3, -P2, -Pi)- 

Note that the order of the spacetime indices in relations ( |7.14j) (as well as in relations 
( |5. 12| ) and ( 5.13Q ) is important. Unfortunately, for the present we cannnot prove similar 



properties for an arbitrary effective subamplitude *T^ S \ since we have not in hand general 
analytic expression for *fX However relations fling ), (ETI3D , and (ff7l4|) for the first 



nontrivial s = 4, 5 provide a reason for assumption that they hold for a general case in a 
form, suggested by Braaten and Pisarski (Eqs. ( |2.3| ) and ( |2.8j ) in Ref. p4[ ) for an arbitrary 
s-gluon HTL-amplitude. 

However as was mentioned at the end of section 5 the Ward identities for effective 
amplitudes have not take place in the form as for usual HTL-amplitudes. The relations 
between higher and lower effective amplitudes can be obtained by direct contraction of 
the expression Q5.6|) and ( [7.13 ) type with momentum variable p. But it is more convenient 



to derive similar relations in approach suggested by Blaizot and Iancu in Ref. || through 
corresponding differentiation of the conservation law for the color current, which in the 
momentum representation has a form 

P%» = -Wf abc j A b ^ Pl )Jfa 2 )5( P - Pl - P2 )d Pl dp 2 . (7.15) 

A gauge field A bfl (pi) on the right-hand side of Eq. (|7.15|) is considered as interacting 
one, i.e. specified by a series (|7.7|) . It is convenient to represent the color current J nip) 
entering into Eq. ( |7.15| ) as a sum of two parts: linear and nonlinear ones with respect to 
interacting field A^p) 

J afl (p) = J?{P) + JTM- (7-16) 

By relation ( |7.6| ) a nonlinear part of the current can be represented in the form of expan- 
sion in terms of a free field 

oo 

JttM(p) - Jn"l[A {0) }( P ) = E J {sW (^°\ . . . , A<°>). (7.17) 

s=2 

By virtue of the fact that a current </^,[^4](p) involves also bare three- and four-gluon 
vertices into its definition, for consistency one need to take a linear part of the current in 
the form 

J?(p) = -*V^(p)A a v (p), (7.18) 

i.e. a coefficient of proportionality between J^ 1 and A afl is here the inverse resumm gluon 
propagator (without gauge fixing term), and not the polarization tensor, as in usual 



28 



definition. Substituting ( |7.16| ) into ( |7.15| ) with regard to ( |7.17p and ( |7.18| ) we result in 
the following expression instead of (|7.15|) 



igf abc A b (p 1 rv^(p 2 )Al(p 2 )5(p-p 1 -p 2 )dp 1 dp 2 



(7.19) 



Wf abc [A b J Pl )J^ L [A^}(p 2 )S(p-p 1 -p 2 )dp 1 dp 2 . 



Here, on the left-hand side we taken into account orthogonality of an inverse propagator 
p^T) -1 ^ (p) = 0. Differentiating Eq. ( [7.19| ) with respect to free field A^ a , considering 
( [7.71 ) and (|7.8|) , one obtains relations between effective amplitudes *T ( - S \ A contraction 
momentum p with the effective amplitude *r( s ) is expressed here in the form of combina- 
tion of the same effective amplitudes of lower order, *r < - s '- ) , s' < s. 



8 Gauge invariance of the effective amplitudes 



The properties ( 5.12 ), (|5.13|) and ( 7.14 ) for effective amplitudes are trivial in the sense 



that they separately hold for every typical group of terms, formed effective amplitudes^]. 
These properties do not reflect more deep connection between different terms involved 
into different groups. Here, we consider one further property of the effective amplitudes, 
for which such a connection completely manifests. 

In proving the gauge-invariance of the nonlinear Landau damping rate, we have shown 
|H that the function ( |5.11|) entering into the definition of the matrix element (|5.7| ) can be 
introduced in its simplest reduced form, 

*f(p, -pi, -p 2 , -p 3 ) = p 2 p\p 2 2 pI{ *r 00 oo(p, -Pi, ~P2, ~Ps) 

- r 00 z,G», -Pi, ~P + Pi) *V vv \p 2 + Pa) *IVoo(pa + Ps, ~P2, ~Ps) (8.1) 
-T 00u (p, -p 3 , -p + p 3 )*V uu ' (p 1 + p 2 )*T u/Q0 (p 1 + p 2 , -p 2 , -pM 

) on— shell 

The proof of validity of this reduction is based on the use of the Ward effective identi- 
ties for HTL-amplitudes and the plasmon mass-shell condition. The right-hand side of 
last expression, as shown in Ref. [Bl is identical both in temporal and covariant gauges. 
The latter is obtained from Eq. (|5.11|) by the following replacements of projector and 
propagator, (see Eqs. ( ^.10[ ) - ( |3.12| )) 



% v {p) -> *V, v (p) = -P, u (p) *A*(p) - Q, u (p) *A\p) + f ZV(p)A°(p), (8-2) 



6 In the case of four-plasmon function (|5.6| ) we look at two groups of terms differ in structure: four- 
gluon HT L-am plitude and group consisting of two terms with three-gluon HTL-amplitude. For five-point 



function (7.13) we have already four groups of terms with essentially different structure. 
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where 



P,, (p) = 9„„ - £V(p) - Qu,{p), QM = "Zl • A °W - . - 



and £ is a gauge parameter in the covariant gauge. All terms on the right-hand side of 
Eq. (|3.1|), which include gauge parameter, vanish on mass-shell. 



We can assume that a similar reduction holds for an arbitrary effective amplitude. 
However a proof of this statement is impossible in the general case by reason of absence 
of general expression for *T^ S \ when a proof by induction is allowed. The only thing that 
we can make is to consider a contraction similar to ( |3.1|) for five-point effective amplitude, 
exact form of which ( [7.13[ ) is known. It should be noted that as distinct from 



given contraction is of no any physical meaning on the plasmon mass-shell, because the 
process of five-plasmon decay is kinematically forbidden. Since here, the intermediate 
calculations are cumbersome, we give only the net result. Using effective Ward identities 
for HTL-amplitudes |fj, |24j] and mass-shell condition, we derive 



T m .../u(P> ~Pi, • • • , ~Pi) ^(p)« Ml (Pi) • • • ^ 4 (P4)| on _ shcll = -P 2 p\p\p\p\ 

T 00 ... Q (p, -pi, -Pi)-{ TWO, -Pl, ~P2, -P3-P4) *^ W '(P3+P4)*f 1 /Q (p3+P4, -P3, ~Pi) 

+ (circular perms, of 2, 3, 4) | 
- \^m v {p,-Px,-p + Pi)*i> vv '{p-PiT^p'mo{p-Pi,-P2,-p^-Pi) (8.3) 
+ T(M)G», ~P + Pi, ~Pi) *V uv '(p - p 4 )*f „>ooo(P ~ Pi, -Pi, -Pi, -Ps)} 
+ ^ou\(p, ~Pi - Pi, ~P3 ~ Pi) *V vu '(pi + p 2 )*f u'mijPi + Pi, ~Pi, -P2) 

X *V XX '(p 3 +p 4 )*f A'00(P3 +Pi, -P3, -Pi)} , „■ 

J on— shell 

As we see, the result is similar to ( |S.lD . It is also easily to check that all terms with a gauge 
parameter in Eq. ( |8.3| ) vanish on mass-shell. If we perform a replacements ( |8.2j ) on the left- 
hand side (i.e. define a matrix element of five-plasmon decay in a covariant gauge), then 
after analogous computations we lead to the same expression on the right-hand side of ( ft.3| ) 
with only distinction in a common signF]. This fact of a gauge non-invariance immediately 
points to the fact that odd effective amplitudes taken on mass-shell of plasma excitations 
are of no physical meaning, as said above. All decay processes involving odd numbers of 
plasmons are kinematically forbidden. Therefore in an expansion of the nonlinear current 
( [7.17 ) all effective currents *jjf' a with even number of free fields (correspondingly 



containing odd effective amplitudes) by 5-functions are equal to zero on mass-shell. 



At the end of this section we note that in spite of the fact that result ( |3.3|) is of only 
pure methodological meaning, nevertheless two examples flB.ip and ( |8.3| ) provide a reason 

7 Such distinct in a sign in calculation of a similar contraction with an effective amplitude with odd 
number of external legs can be already found for first three-point effective amplitude T W1AI2 [p, —pi, —P2)- 
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to use of considerably simple expressions of ( |8.1| ) type for all (2n + 2)-matrix elements 
jypai...a 2n+ i j n particular calculations. 



9 The Vlasov-Boltzmann equation for color plasmons 

In a previous sections we have considered the problem of deriving kinetic equation of a 
Boltzmann type, describing decay processes involving colorless plasmons. The number 
density of plasmons here has a trivial color structure: N p ab = 5 ab N p . Such a structure of 
the plasmon number density takes place, when there is no external color current or/and 
color mean field in the system. Now we assume that a time-space dependent external per- 
turbation (e.g. external color current j^ ta (x)) starts acting on the system. In the pres- 
ence of the external color perturbation, the soft gauge field develops an expectation value 
(v4^(x)) = A"(x) 7^ 0, and the number density of the plasmons acquires a non-diagonal 
color structure. Thereby, we lead to the problem of construction of more complicated 
kinetic theory for color plasmons. The kinetic (matrix) equation in this theory will take 
into account such purely non-Abelian effect as the precession of the color charge of the 
plasmon, the existence of which qualitatively differs non-Abelian plasma from Abelian 
one, where the plasmons do not carry electric charge. Our further consideration will 
be to a certain extent of phenomenological character, being only a first step towards to 
construction of a kinetic theory for color plasmons in hot QCD plasma. For rigorous 
justification of all assumption it should be invoked the methods of an off-equilibrium field 
theory. 

The assumption that a structure of the kinetic equation for color plasmons is similar 
to the structure of the kinetic equation for hard transverse gluons in the form proposed by 
Arnold, Son and Yaffe in Ref. ||25|| , will be our initial theoretical premise to construction 
of above-mentioned transport theory. More precisely, we expect the time-space evolution 
of N p = (N p ab ) to be described by 

(V t + V l p ■ V x )Nl - \{(E{x) + (V p x B(x)))., V P X} = - C [N l p ], (9.1) 

where is a covariant derivative acting as 

V^Nl^d^ + igiA^x),^}, 

A^ = T a A° il is mean soft-gluon field expressed in terms of Hermitian generators in ad- 
joint representation T a ((T a ) bc = -if abc , tr(T a T b ) = N c 5 ab ) } with [,] and {,} denoting 
the commutator and anticommutator in color space, respectively; £ l {x) and B l (x) are 
the mean chromoelectric and chromomagnetic fields. Futhermore we consider that the 
collision term C [N p ] has a following structure 

C [N l p ] = \ {N l p , T d [iV p ]} - \ {(1 + N l p ), r ; [iV p ]} - . . . , (9.2) 
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where TdfiVp] = (T^'fiVp]) and ri[iV p ] = (r? a ' [N^]) represent a generalized decay and re- 
generating rates of color plasmons, respectively. These rates can be formally represented, 
like a colorless case in the form of functional expansion in powers of the number density 
of color plasmons 

oo oo 

rrK] = E rf +1)aa K], rfWl) = £ rp^-'ijvj], (9.3) 

n=l n=l 

where 

r^ aa '[Nl] = J W^(P, Pl, .. , Pn; Pn + 1, • • • , P 2n+l) <^ ■ ■ • JV#"» 

(9.4) 

X (1 + A^ +i )<^< +1 . . . (1 + < n+ J^^4n +1) 

r [*+i) = J dT ^w^(p, Pl , . . . , Pn ; Pn+1 , . . . , p 2n+1 ) (1+iV^rW . .(l+JVjJ*» 

xiViS 1<4 - 1 ...JViK 1,C4 - 1 - (9-5) 



Here, w^ 2 } = ^"X-^+iS^i-^+^p, Pl , . . . , Pn ; Pn+1 , . . . , p 2n+1 ) is the scattering prob- 
ability defined as 

{a'; a}/ \ 
W 2n+2 IP) Pl) • - - ) Pn! Pn+1) - - - , P2n+lJ — 

= M* {a ' } (p, P i, . . . , p n , -Pn+i, • • • , - P 2„+l)M {a} ( P , P i, . . . , Pn , - P „+l, • • • , -P2n+l) 

(PiiOij) (Pi 1 a< 1 :Pi 2 a <2) (piai,...,PnOn) 

+ ^ M* {a ' } M {a} + ^ M* {a ' } M {a} + . . . + Y, M* {a ' } M {a} 

l<il<n l<ii<i2<n 

with the same matrix elements M\ a * (p, P i, . . . Pn , — Pn+ i, . . . , — P 2n+i) as was defined by us 
in colorless case by general relation ( |4.15| ) (see, however the end of this section). Here, the 

summing symbol Z)i<ii<n designates that besides summing over all possible interchange 
of momenta p i± <-> —Pj x , where indicies z 1; ji run 1 < i% < n, n + 1 < ji < 2n + 1, it needs 
to interchange and color indices a^, «-> at one time. The other summing 

symbols are taken in a similar fashion. In arrangement of color indices in Eqs. ( |9.4|) 



and ( |9.5| ) we follow rule proposed by Arnold, Son and Yaffe in Ref. [53] in the context 
of deriving collision term for hard color particles. The dots on the righ-hand side of 
Eq. ( |9.2| ) is refered to contributions containing commutators of [Ren fl , N^] type, where 
n R is retarded soft-gluon self-energy. 

In the case of four-plasmon decay, the scattering probability w\ a ' a ^ by general expres- 
sion has a form 

{a'; a} 



< ;a| (p, Pi; P2, Ps) = M* a '< a '^(p, Pi, -P2, -p 3 )M™ 3 (p, Pl , -p 2 , -p 3 ) 
+ M* a ' a '^ ( P , - P2 , Pl , - P3 )M™ 3 (p, - P2 , Pl , -p 3 ) 

+ M* a ' a '^<(p, -p 3 , - P2 , P l)M™ ai (P, - P3 , - P2 , Pl). 
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By direct calculation, using an explicit expression ( |5.5| ) for color structure of an effective 
amplitude *Y\ a ^^ and properties for subamplitude *f\^ (|5.12| ), ( f5.13| ), it can be shown 
that a sum of two last terms in this expression is equal to doubled value of a first term. 

It needs to supplement the Vlasov-Boltzmann equation ( |9.1|) by mean field equation, 
defining a change of mean field in a system in self-consistent manner 



V^x)^(x)=jf^(x)+f^i 



X) 



where induced current j£ lasm (x) = (j£ lasm (x), j plasm (x)), 

Jo Plasn » = 9T a j tr(T a N l p ), f^(x) = gT a J M tl (T a N l p ) (9.6) 

is the color current resulting from color-plasmon number density and j^ xt is the external 
current plays a part of initial color perturbation. 

It is clear that dynamics of nonlinear interaction of soft color excitations is considerably 
more intricate as compared with colorless ones. To form a certain notion of this fact we 
consider in detail the linearized version of Vlasov-Boltzmann equation keeping only 

the first term, n = 1, in expansions (|9.3|) . We write the number density of color plasmons 

as 

N p ab = N l eq (p)S ab + 5N p ab , 

and parametrizate the off-equilibrium fluctuation of the occupation number similar to 
Eq. ( |5.8| ), where now the function W l p is a color matrix in the adjoint representation: 
Wp = W p a T a . After some cumbersome algebraic transformations, we derive a linearized 
kinetic equation for the color plasmons from Eqs. (9.1) - (9.5) 



(V t + V p .V x )W p = -g (V p • £{x)) - 5C [W l p ], (9.7) 
where linearized collision term 5C has a following form 

A - P hA ,ii f rf Pi d Pz d p 3 4 , , , , 



dpi dp 2 dp 3 

3(2tt) 

<(pi)«( P2 ) + i)«(p 3 ) + i; 



X *CP + P. - P 2 - P3) — TOTI) 



X 



+ Re(T 1 T 2 l{K- l(K,+K 3 )})- 



2 

On the right-hand side of Eq. ( |9.8|) we introduce for brevity a notation 
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Here, the function *T is defined by Eq. ( |5.11| ), and a function T 2 is defined from T\ by a 
replacement 

*f(p,pi, -p 3 , -P2) -> *f(p, -P2,Pi, -Pa). 
In deriving ( |9.8| ) we use the following identities: 

tr{T a T b T c ) = ~ iV c / aftc , (T a T b ) ti{T a T c T b ) = - N 2 T C , (T a T b T c T a T b ) = 0, etc. 

Since the equilibrium plasmon number density is proportional to the identity, the com- 
mutator terms on the righ-hand side of Eq. Q9.2p vanish. Therefore if we assume that the 
off-equilibrium function SN l is perturbative small, then a linearized Boltzmann equation 
in a certain sense is exact. 



The linearized collision term (|9.8|) for color plasmons has noteworthy structure (cp. 
collision term for colorless case (|5.9Q , (|5.10|) ). The expressions in curly brackets behind 
squared moduli |TjJ 2 and \T 2 \ 2 coincide in a structure with a similar expression entering 
into a linearized collision term for color fluctuations of number density of hard gluons in 
form which was proposed by Blaizot and Iancu jnj] . Besides, here additional " interference" 
term Re (TiT 2 ) appears. Going immediately after an expression in curly braces does 
not coincide with none of known expressions. In spite of some complexity of obtained 
expression for 5C [Wi] it provides remark internal symmetry, directly connected with 
existence of relation (|5.12| ). Let us rewrite this relation in terms of T functions ( |9.9|) . For 
this purpose we introduce a third function T 3 , that is defined from 7{, by replacement 

*f (p,Pi,-p 3 ,-P2) -> *f (p,px, -p 2 , -Pa)- 
Then from (p-12[) it follows 



r 1 + T 2 + % = 0. 

Now we set, for example, a function T 2 = —T\ — % and substitute it into expression in 
round brackets in integrand (|9.8|) . After simple algebraic transformations we have 



|7"iI 2 { w p- iKi-liKi + K^ + w^K- iKs-lo^k+KJ} ( 9 - 10 ) 

The structure of the last expression with an accuracy of replacement of momenta fully 
coincides with a structure of Eq. ( |9.8| ). The linearized collision integral Q9.SQ may be 
rewritten in terms of functions T 2 and T3 by setting T\ = —T 2 — T3. We obtain similar 
expression with another permutation of momenta. Adding such obtained three expressions 
and dividing by 3, we result in an expression that is more symmetric relative to the external 
soft momenta and will be used somewhat below 

Q(p,Pi;p 2 ,p 3 ) = \{\T 1 \ 2 {2W l p - \w l v:i - \{W^ + w l P2 )} 
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+ \T 2 \ 2 {2W l p - \ W l pi - ~ (W< 2 + W< 3 )} + \%\ 2 {2W l p -\K,-\ (Ki + KJ} 
+ Re (T^*) {W l p - \ (W l pi + W l P3 )} + Re [%T*) {w l p - \ (W l p2 + W l j} (9.11] 
+ Re(T 2 T 3 *){w< - Uw l pi +W l p2 )}). 



By direct calculations we show that the color current induced by off-equilibrium fluc- 
tuations of the color plasmons number density 5N p ab , covariantly conserves. In terms of 
functions W l p a , the expressions for components of color current ( |9.6| ) is rewritten in the 
form 

Plasm, ^__ nT af ^ f dN oM\, A ,la jpla-m M _ _ f dp fdN^jp)}^ 



(9.12) 



By linearized Boltzmann equation (|9.7|) we have 



(9.13) 

xN l eq (p)N l eq ( Pl )(l + <( P2 ))(1 + i\^(p 3 ))Q(p, Pl ; p 2 , Ps ), 

where the function Q on the righ-hand side is defined by Eq. ( |9.11| ). Now we show that 
integral on the righ-hand side of the last expression equals to zero. As a first step we 
symmetrize integrand in Eq. ( |9.13| ) with respect to permutation p Pl . Using definitions 
of functions T ly T 2 and T 3 and property of invariance of the function *T(p, —pi, — p 2 , —pz) 
(Eq. (|5.13| )), when the momenta order is reversed, we obtain laws of transformation of 
these functions with replacement p <-» pi 

Ti^Ta, T 2 ^T 2 , (9.14) 

Furthermore, we replace p pi in expression (" [9 .11 ) with regard to ( 9.14Q . Adding (|9.11|) 
with such an obtained expression and dividing by 2, we lead to the expression symmetric 
over replacement p «-> pi: 



X 



As we see, this expression is automatically symmetric and with respect to permutation 
P2 ^ P3 (here, transformation laws of functions T coincide with (]9.14 )). 

Let us consider now a crossed symmetry p <-> p 2 ,pi <-» P3. The statistical factor in 
integrand ( |9. 13| ) is symmetric by virtue of identity 

<(p)<(Pi)(l + <(P 2 ))(1 + <(p 3 )) = (1 + N l eq (p))(l + <(pi)X(p 2 X(p 3 ). 



35 



The transformation laws of functions % are trivial in this case 



By using these transformation laws and expression ( |9.15| ) we see that integrand on the 



right-hand side of Eq. ( |9.13| ) is odd for replacements p <-> P2, pi <-> P3. By virtue of this 



fact the integral equals zero and a color current (|9.12|) is covariantly conserves. 



Closing this section note that without rigorous transport theory for color plasmons in 
hand here, we given no any estimations for the amplitude of mean gauge field A^. We 
can only made proposal that for sufficiently strong background field, a structure of the 
kinetic equation (|9 . 1| ) can be more involved. The interaction matrix elements y[ aai " M2n+1 
(or effective amplitudes), for example, can depend by itself in a highly nontrivial way on 
mean field. In other words, they can contain insertions of external mean field, in general, 
of arbitrary orderQ. This leads to the fact that mean gluon field will be "involved" in the 
dynamics of color plasmons not only through the left-hand side of the kinetic equation, 
but through a collision term. 



10 Conclusion 

On the basis of pure gauge sector of the Blaizot-Iancu equations we have obtained the 
transport equation for the colorless plasmons, taking into account (2n + 2)-plasmon de- 
cays. The algorithm of successive calculation of the probabilities of the decay processes is 
proposed. The limiting value of the plasmon occupation number iV p (~ 1/g 2 ) is defined, 
whereby the decay processes with an arbitrary number of soft external legs give contribu- 
tions of the same order in the coupling constant to the change of iV p , and transport theory 
becomes essentially nonlinear. Here, it should be pointed to certain difficulty of a prin- 
cipal character, which arises when we going over from low excited state (A^X) ~ \/gT) 
to highly excited one (A^X) ~ T) of system. 

The kinetic equation (|2.2|) describes the processes of plasmon scattering. However 
moreover, there is a collective plasmon interaction resulting in frequency shift of plasmons 

^ P -> ^ P = ^ P + A^ p . 

Here, u l is the frequency renormalized due to the interaction. In the linear plasma theory, 
i.e. in the theory with infinitesimal amplitudes of plasma oscillations, the frequency 
spectrum represents the function of values (T, g, . . .), not dependent on waves energy (or 
their occupation number). In the framework of nonlinear theory, the nonlinear dispersion 
correction represents in general case a functional of plasmon number density iV p . To 

The usual HTL-amplitudes in the presence of a background gauge field, proposed by Blaizot and 
Iancu in Ref. H can be an example of such a dependence. 
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lowest nontrivial order in iV*, the nonlinear frequency shift of longitudinal oscillations is 
defined by the formula 



p y C J (2vr) 3 Pl V 2u l p J\ 2u l pi 



u 2 (p)u 2 (pi 



-Re*f(p,px, -p, -pi) 



on— shell 



+ o((N p y) + ..., (10.1; 



where the function *T(p,pi,—p,—px) is defined by Eq. ( |5.11|) . At the soft momentum 



scale, for the thermal fluctuations of a gauge field we have an estimation for the correction 

Au l p ~ g 2 N c T < J p . 

Such for low excited state the frequency shift is perturbative small. In a limiting case of 
a highly excited state (6.4) from Eq. (|10.1|) we have the following estimation 

Au l p ~ gN c T ~ u l p , 

i.e. the nonlinear frequency shift becomes of the same order[] as a linear part of a spectrum 
uj l p . The spectrum broadening over frequency takes place, and for strong field < T 
a dependence I l (p,x) on a frequency (which in this work is defined in the form of the 
quasiparticle approximation ( |4.12|) ) has no connection with 5-function. Therefore, the 
kinetic equation in the form ( |2.2j ) in this limiting case becomes in general inappropriate 
for description of nonlinear dynamics of soft excitations, and here, it should be derived 
more general equation (or system of equations) for initial Wigner function I Hp, x), taking 
into account the spectrum broadening effects. The examples of conctraction of such an 
equations for usual plasma can be found in Ref. |2"B| . Nevertheless, we suppose that 



somewhat simplified approach to kinetic description of highly excited state of hot gluon 
plasma considered in this paper will be useful in at least qualitative analysis of complicated 
dynamics of interaction of soft boson excitations taking place in the extremal medium 
state. 

At the end of this section we note that at present we study a posibility of an alternative 
way of description of nonlinear plasmons dynamics based on a Hamiltonian formalism. 
The cornerstone of this approach is a fundamental fact that the equations describing a 
collisionless hot gluon plasma in HTL-approximation, possess a Hamiltonian structure, 
which was proposed by Nair, Blaizot and Iancu |7|, |8|, |2"Tfl . The last circumstance enables 
us to developed (at least for low excited states) independent method of derivation of 
kinetic equation for soft gluon plasma excitations. Within the framework of Hamiltonian 
approach, matrix elements of (2n + 2)-plasmon decays are obtained as result of a special 
canonical transformation, simplifying the interaction part of the plasmon Hamiltonian. 
The exact consideration of this approach will be subject of a separate publication. 



3 In this case to be sure in the right-hand side of Eq. ( |l0.l| ) all terms of higher power in are of the 



same order in g, as the first term (similar to expansion of generalized rates (2.3)) 
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FIGURES 



FIG. 1. The matrix element for four-plasmon decay. The wave lines denote soft quasiparti- 
cles (plasmons) and the blob stands for HTL resummation. 

FIG. 2. The typical tree-level Feynman diagram for (2n+2)-plasmon decay with amputate 
external legs. 
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FIG. 2. 



